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Critical properties in the dynamical charge correlation 
function for the one-dimensional Mott insulator are studied. 
By properly taking into account the final-state interaction be- 
tween the charge and spin degrees of freedom, we find that 
the edge singularity in the charge correlation function is gov- 
erned by massless spinon excitations, although it is naively 
expected that spinons do not directly contribute to the charge 
excitation over the Hubbard gap. We obtain the momentum- 
dependent anomalous critical exponent by applying the finite- 
size scaling analysis to the Bethe ansatz solution of the half- 
filled Hubbard model. 

PACS: 75.10.Jm, 05.30.-d, 03.65.Sq 



I. INTRODUCTION 

Recently strongly correlated electron systems in one 
dimension (ID) have attracted much interest. In par- 
ticular, charge and spin excitations for the Mott insu- 
lator have been intensively studied. Various experimen- 
tal methods, by which the dynamical correlation func- 
tions are directly observed, have revealed striking prop- 
erties in the Mott insulator. For example, the recent 
photoemission experiments [^-^| for the ID compounds 
SrCuC>2, Sr2CuC>3 and NaX^Os have clarified dynamical 
properties of a single hole doped in the ID Mott insula- 
tor. These experiments have indeed stimulated extensive 
theoretical studies on the one-particle Green function for 
ID correlated systems, which have been done numeri- 
cally §^ and analytically. fl|fH19| The overall feature 
of the spin-charge separation found in the photoemission 
experiments has been explained rather well. jq| In this 
connection, it has been claimed that the final-state inter- 
action between the charge and spin degrees of freedom 
plays the crucial role to determine the edge singularity in 
the spectrum, based on the conformal field theory (CFT) 
analysis. 

The dynamical charge correlation function is another 
key quantity to explore the characteristics of the Mott 
insulator, which is particularly important to analyze op- 
tical experiments. This quantity has also been stud- 
ied intensively so far for ID correlated electron systems. 
p(i|-p5| For example, the recent theoretical treatment of 
a metallic system close to the Mott transition [E2j has 
clarified that the change of the weight in the dynamical 
charge correlation function clearly describes the charac- 
teristics of the metal-insulator transition. 



Motivated by the above investigations, in this paper 
we study critical properties in the dynamical charge cor- 
relation function for the ID Mott insulator. For this pur- 
pose, we consider the ID half-filled Hubbard model as a 
Mott insulator and use the exact solution of the Bethe 
ansatz method. pfj|-|28|] Applying CFT techniques, we an- 
alyze the exact finite-size spectrum to obtain the critical 
exponents for the dynamical charge correlation function 
exactly. We clarify how the edge singularity in the mas- 
sive charge excitation spectrum over the Hubbard gap 
is controlled by massless spinon excitations by properly 
taking into account the final-state interaction between 
the charge and spin degrees of freedom. In particular, we 
point out that a particle-hole charge excitation is consid- 
ered to act as two mobile impurities in massless spinons, 
and the resulting scattering phase shifts are the essential 
quantities to determine the edge singularity. 

This paper is organized as follows. In § 2, based on the 
exact solution of the ID Hubbard model, [^6| we present 
the basic formulation following refs. 27 and 28 to inves- 
tigate a particle-hole excitation over the Hubbard gap. 
Then in § 3 we calculate the exact finite-size spectrum, 
and derive the edge singularity in the charge correlation 
function by employing the finite-size scaling idea in CFT. 
We discuss the anomalous critical behavior in the sin- 
gularity by evaluating the momentum-dependent critical 
exponents. Brief summary is given in § 4. 



II. CHARGE EXCITATIONS IN THE MOTT 
INSULATOR 

Let us start with the ordinary ID repulsive Hubbard 
model, 
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where cj is the creation operator for electrons and h is 
a magnetic field. We henceforth set t = 1 for simplicity. 
The exact solution of this model |2(| is given by a set of 
the Bethe equations for the charge rapidities kj and the 
spin rapidities A Q , 

M 



Nkj = 2ttJj - ^2 2tan~ 1 [ 

0=1 



sin k-i — A 



0} 
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jr 2tan- 1 [ ^ Zph ] = 2 -Kj a + 2tan- 1 [^^], 

where we set c = [7/4 (> 0). jV e and M represent the 
number of electrons and down spins, and N is the number 
of the lattice sites. For the repulsive case U > 0, the 
rapidities are real numbers for the ground state as well as 
for excited states within the lower Hubbard band, which 
are classified by the set of the quantum numbers Ij and 
J a for the charge and spin sectors. 

A. Particle-hole excitations at half-filling 



which are recast into more convenient form, 

A = -(sinfe; + smk m ). (5) 

Note that the unknown parameters for excitations are 
now reduced to two rapidities ki and k m thanks to the 
above constraints, which naturally reproduces the fact 
that a particle-hole excitation is classified by two inde- 
pendent momentums. Eqs. (|^) and (|5|) are our start- 
ing equations to analyze the finite-size corrections to a 
particle-hole excitation over the Hubbard gap and to get 
further information on the charge correlation function. 



For the half-filled case (N e = N) the available quantum 
numbers Tj for the charge sector are completely filled, 
giving rise to the incompressible state of the Mott in- 
sulator. So, charge excitations over the Hubbard gap, 
which are necessary for the evaluation of the dynamical 
charge correlation function, cannot be described by real 
rapidities, 

It is known that the charge excitations of particle-hole 
type are specified by complex charge rapidities in bound 
pairs. |2?J Following the methods of Woynarovich, psfl 
we briefly summarize how to treat such complex rapidi- 
ties. As a simple charge excitation at half-filling, we here 
introduce a pair of complex charge rapidities, which are 
denoted as k^ — k ± i\ with the real quantities k and x- 
In the thermodynamic limit, these charge rapidities are 
coupled with the corresponding spin rapidity A to yield 
a bound pair, |27| 



sin(ft ± ix) — A T ic 



(2) 



Since this pair adds two particles to the system, we need 
to introduce two holes, labeled by the real charge ra- 
pidities ki and k m , in order to keep the particle number 
unchanged. We thus obtain the Bethe equations for a 
particle-hole excitation over the Hubbard gap, p8| 
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B. Dispersion relations in magnetic fields 

Before computing the finite-size corrections, we first 
summarize the results on the dispersion relation by 
slightly extending the method of ref. 28 to include the 
effect of magnetic fields. Following a standard way, let 
us introduce the distribution functions for the charge and 
spin sectors, which are denoted as p(k) and c(A), respec- 
tively. 

1 dA 

p(k) = VI — cos kK\2 (sin k — A)ct(A) 

27r J-x 27r 



1 



2ttN 

— K 2 i(X- sin k)p(k) - -^—^ ^ #2i(A-sinfc Q ) 



cos kK\ 2 (sin k — A) 
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dA 



o-{\)+ I —K 22 (\~\)a(\ 



(6) 
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where K 12 (x)=K 2 i(x) = 2c/{c 2 + x 2 ), K 22 {x)=Ac / (Ac 2 + 
x 2 ). Note that a particle-hole excitation is incorporated 
in the integral equations through the terms of order of 
1/N including ki and k m . 

The energy and momentum for this excited state are 
given by, 

M 

E = V" (h c - 2 cos kj) + h s 

j^l,m a— 1 

+ (h c — 2cos/c + +h c — 2cosfc~), 



These equations are subject to the constraints for the 
momentum conservation, 



p= k j + ( k+ + k ~)< 



where h c = —h/2 and h s — h. Using (^), we recast them 
to simple formulas, 



oj c = E - E = U - e c {h) - e c (k m ), 
q c = P - Pa = -p(k) - p(k m ), 



(8) 
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where Eq and Pq is the energy and momentum of the 
ground state. Here we have introduced the dressed en- 
ergy e c {k) and the dressed momentum p(k) for the sim- 
plest charge excitation (often referred to as holon), which 
are given by 



£c(k) = — 2 cos k — 2 / Ji(w) cos(w sin k)duj 

y ' J Q ujcoshcuj v ; v ; 

1 7T 

■cosh -1 , , . — ; — e s (A)dA, 



p(k) = k + 



a|>a 4c 2c(sinfc - A) 

^ — CUJ 

-Jq(uj) sin(w sin fc)dw 



i cosh c 



A|>A 



G(sinfc - A)ct(A), 



(9) 



where J n {iio) is the n-th order Bessel function. This 
charge excitation is coupled with spin excitation 
(spinon) , whose energy e s is given by the solution to the 
following integral equation, 



/■Ao j\' 

e s (\)=s°sW- / -^K 22 (X-X)e s (X) 



2tt 



(10) 



where, 



G(A) 



h-2 



dk 
2^ 
sin uj\ 



' cosh < 



cos 2 kKi2 (sin k — A) , 



■duj. 
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The dispersion for the particle-hole excitation spectrum 
is described by w c (g), which includes the effects of mag- 
netic fields. 

In Fig. 1 we show the charge excitation spectrum (q = 
q c ) for a given U in several choices of magnetic fields. The 
charge excitation spectrum distributes in the continuum 
over the Hubbard gap. It is seen from the figure that the 
lower edge of the spectrum changes its character at the 
critical momentum q — it ± q L (h). For n — q L < q < 
Tt + qL, the lower edge is given by changing a rapidity k m 
with the other being fixed as ki = —tt. This excitation 
corresponds to a particle-hole excitation from the top 
of the lower Hubbard band to an excited state in the 
upper Hubbard band. On the other hand, for < q < 
n — qL and tt + q^ < q < 2tt, the lower edge is featured 
by the excitation with ki — k m , which makes a particle 
(hole) with the momentum q/2 (—q/2) in the the upper 
(lower) Hubbard band. So, the corresponding critical 
behavior of the edge singularity is different between these 
two regimes, which will be explicitly studied in the next 
section. 



FIG. 1. Charge excitation spectrum for U = 6 in several 
choices of magnetic fields: (a), (b) and (c) correspond to the 
case of the magnetization m = 0,0.5 and 1. The excitations 
are allowed for the hatched region. The lower edge of spec- 
trum changes its character at the critical momentum n ± qL- 
In the case of U — + oo or m — > 1, qL equals to in < q < 2tv. 
On the other hand, qL becomes n at U = 0. 



III. CRITICAL BEHAVIOR IN DYNAMICAL 
CHARGE CORRELATION FUNCTION 

We now consider the dynamical charge correlation 
function for the Mott insulator and show that the anoma- 
lous power-law behavior appears in the spectrum, which 
is controlled by massless spinon excitations subject to the 
two phase shifts due to a particle-hole excitation. 

Let us start with the charge correlation function, 



D(x,t) = < n(x,t),n(0,0) > 

= ^|<0K0)| M >| 2 c 



iAP u x-iAE u t 



e^-^DsOM), 

{h,k m } 



(12) 



where n — ~^2 a n a is the charge density operator. In the 
third line of the equation, we have approximately sub- 
stituted ui c (ki,k m ) and q c (ki,k m ) given by (||) for the 
energy and the momentum of excited states. Since the 
action of the charge density operator n creates not only 
particle-hole charge excitations, but also induces low- 
energy spin excitations, the contribution from the spin 
sector to < n(x, f)n(0, 0) > is denoted by the spin corre- 
lator D s (x,t). In the following we will mainly deal with 
the Fourier transform of the correlation function, 

D{q,uj) ^Y,{ k ,,k m }J dx I dt e~ i( - q ~ qc)x+Kw ~ We)t D s (x,t), 



for which the main spectrum is featured by the massive 
charge excitation tu c (q) while its critical behavior is es- 
sentially determined by the spin correlator D s (x, t). It is 
now clear that the long-time or large- distance behavior 
of D s (x,t) is important to discuss the critical behavior 
of the dynamical charge correlation function. Therefore 
our remaining task is to study low-energy properties of 
D s (x,t) exactly when a particle-hole excitation is cre- 
ated. 



A. Finite-size spectrum and conformal properties 

In order to apply the methods developed in CFT |29] to 
the evaluation of D s (x, t), let us first compute the finite- 
size corrections to the excitation energy. |3^^l[l We recall 
here that in the calculation of the spectrum in the previ- 
ous section, we have not taken into account the fact that 



(13) 
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massless spinons are scattered by two dressed particles 
created for the charge part (holons) , and have discarded 
the resulting scattering phase shifts in the thermody- 
namic limit. However, it turns out that this coupling be- 
tween the charge and spin degrees of freedom is essential 
to determine anomalous critical properties in the charge 
correlation function. Since we are now considering the 
situation that a particle-hole pair excitation is suddenly 
created, this interaction may be regarded as the so called 
final-state interaction. We show explicitly the above fact 
by correctly evaluating the final-state interaction. 

We can perform the calculations of the finite-size spec- 
trum by extending those previously done for the photoc- 
mission spectrum. fl9f| In order to deal with the finite- 
size effects, let us rewrite the coupled integral equations 
of distribution functions (|^), 

1 dA 

p{k) = h / — coskK 12 (smk - A)ct(A) 

2tt J x - 2tt 

1 



which is referred to as the surface energy in boundary 
CFT. Now, according to the finite-size scaling idea in 
CFT, the scaling dimension x for the spin sector 

is read from the 1/N corrections to the excitation energy, 
which is obtained as, 



1 

4g 



■ n c {k m )f + g(d c (ki) + d c (k m )Y 



(17) 



This formula for the scaling dimension is typical for 
Tomonaga-Luttingcr liquids classified by c = 1 CFT. 
Here, the quantity £ s = £ s (Ao) (referred to as the dressed 
charge) , 



6(A) 



1 



v ° dA 



— X 2 2(A-A)6(A), 



(18) 



2ttA 

— K 2 i(\ - sinfc)p(fc) - — — V" K 2 i(X-wak a 
_ n Ztt ZirlM 



cos kKi2 (sin k — A) , 



features the U(l) critical line of c = 1 CFT when we 
change the strength of interaction or the magnetic field. 
A remarkable point in (|l7]) is that there exist two kinds 
of phase shifts n c and d ci which are caused by the final- 
state interaction between the charge and spin degrees of 
freedom. These phase shifts are explicitly obtained as 



A" 1 



dA 



cr(X)+ I —K 22 (X-X)a(X), 



( 14 ) d c (k a ) 



<7 a (A), 



where we have introduced the asymmetric cutoffs for spin 
rapidities, X^, which are essential to include the final- 
state interaction. Namely, two massive holons created 
give rise to the phase shifts in massless spinons, which 
naturally induce the change in cutoffs in the order of 1/N 
(compare them with (0)). The importance of this effect 
has been previously noticed for the photoemission spectra 
in the Mott insulator. |lj|[ll| In this way, although it is 
naively expected that the massless spinons do not directly 
coupled with the charge excitation, the edge singularity 
in the charge correlation function is affected by spinons 
via the final-state interaction. 

Let us now exactly analyze the excitation energy, in- 
cluding the contribution from low-energy massless spinon 
excitations. By applying the computational techniques 
presented in [jl9| to eqs.(p^), we straightforwardly end 
up with the following formulas, 




°~ a (A) 



<7«(A) 



(19) 



where 



CT a(A)=— AT 2 i(A - sinfc Q ) 
Ztt 



v ° dA 



A 



2tt 



K 22 (X-X)a a (X). 



(20) 



AE = Ae c + Ae s 



= U- e c {h) - e c (k m ) + ^(x + N + + N_), 



where v s is the velocity of massless spinons, 



el (A) 
2na{X) 



A = An 



(15) 



(16) 



It is thus seen that although the scaling dimension x 
for massless spinons is typical for c = 1 CFT, |32 33 
the massive charge sector also contributes to x via the 
phase shifts n c and d c . In this sense, (0) is classified 
as shifted c = 1 CFT, whose fixed point is different even 
from that of the static impurity problem, as pointed out 
by Sorella and Parola. |l6| Since in the present case, these 
phase shifts are considered to be caused by the two holons 
created, the fixed point belongs to Tomonaga-Luttingcr 
liquid with two mobile impurities. p4[ It will be shown 
that these phase shifts are the key quantities to control 
the anomalous low-energy properties in the correlation 
function. 

Exploiting finite-size scaling techniques in CFT, we can 
now write down the asymptotic form for the spin corre- 
lation function as 



D s {x,t) 



1 



(21) 



Note that the first term of the order of unity in ( |15| ) rep- 
resents the energy u c for the massive charge excitation, 



(x- v s t) 2A °(x + v s t) 2A " 

where are conformal dimensions for the spin sec- 
tor, which are related to the scaling dimension as x = 
+ A~ . Thus we can read the low-energy behavior in 
D s (x,t) correctly from the scaling dimension (p"7j). 
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B. Anomalous critical behavior 

Let us now discuss the critical properties in the dy- 
namical correlation function by substituting the above 
results for D 8 (x, t) to (|l3|). As noted in the previous sec- 
tion, there appear two different lower edges in the charge 
spectrum depending on the momentum q, so we deal with 
two cases separately. 

1. Edge singularity for ty — qi, < q < ty + qi, 

We now discuss the critical properties along the lower 
edge of particle-hole excitation spectrum. Let us start 
with the case for ty — q L < q < ty + qL shown in Fig. 2. 
As shown in the previous section, the lower edge is given 
by u) c (q) = U — e(— ty) —e(k m ). In this case, if we neglect 
the contribution from the spin sector, there does not exist 
the power-law edge singularity in the charge correlation 
function, so that the anomalous behavior around the edge 
may be completely determined by spinon excitations. So 
the Fourier transform of the charge correlation function 
is given by 

D(q,u>) oc 6(u - LJ c (q))(tJ - io c {q)) x ^\ (22) 

for oj ~ tL> c (q). The corresponding critical exponent 
Xi (q) directly follows from the scaling dimension for the 
spin sector as 

X 1 (q) = 2x-1, (23) 

where x is given by ( |l7| ) . We emphasize here that the crit- 
ical exponent is dependent on the momentum q. Recall 
that this anomalous behavior for the critical exponent is 
due to the two phase shifts in x, which is caused by the 
final-state interaction. In Fig. 2, we show the results for 
X\ (q) as a function of q for a given magnetic field. In the 
limit of vanishing magnetic fields, the exponent gradually 
approaches the value of —1 because the scaling dimension 
x tends to zero. On the other hand, exactly at ft- = 0, 
the power-law singularity is expected to disappear, which 
implies that X\ — 0. We can show that there indeed 
exists an energy scale which characterizes the crossover 
between these two behaviors. Namely, only in the small 
region close to the edge, (u — uj c (q))^ 2x ^ 1 " > < 1/T(2x), 
we can observe a power-law singularity with the expo- 
nent X\ — 2x — 1 where T is the gamma function. When 
the magnetic field is decreased, this region is gradually 
shrunk because r(2ar) — > oo, and is eventually replaced 
by D(q,u> — > u c ) ^constant. 

2. Edge singularity for < q < ty — qL and ty + qL < q < 2ty 

Let us now turn to the momentum region < q < 
ty — qL and ty + qL < q < 2tt in Fig. 1, in which the lower 



edge is given by the excitation with ki = k m . Along 
this edge there occurs a specific excitation where a par- 
ticle (a hole) is created at symmetric points in the upper 
(lower) Hubbard band with the momentum q/2 (— q/2) 
for a given q. So, we should be a little bit careful to 
deduce its critical behavior. We start by discarding the 
spinon contributions to clearly see the origin of the sin- 
gularity. Let us first recall that for a given q the charge 
excitation energy lu c in (^) takes its minimum at ki = k m 
when we change the values of ki and k m . So, around the 
edge, lo c has the following properties for a fixed momen- 
tum q, 

uj c \ km =h = 0, uj c \ km=kl > 0, 

(24) 

which naturally leads to the well known square-root be- 
havior like the van Hove singularity, [^2j (ui — ujc(q))^ 1 / 2 . 
This singularity depends on neither the coulomb inter- 
action U nor the magnetic field h, so long as the spinon 
contribution is neglected. 

Now let us take into account the final-state interaction 
between the charge and spin degrees of freedom. Then 
low-energy spinons subject to the phase shifts show up in 
the infrared singularity. By properly incorporating this 
contribution, we finally end up with the edge singularity 
for < q < ty — qr, and tt + q^ < q < 2tt, 

D(q,Lo)Yx9(u-u c (q))(u-u c (q)) x ^\ (25) 

where the corresponding critical exponent is 

X 2 (q) = -± + ±(2x-l). (26) 

The first term of —1/2 comes from the square- root singu- 
larity in the dispersion relation, while the latter contri- 
bution reflects low-energy spinon excitations. In Fig. 2 
we show the momentum dependent critical exponent for 
< q < ty — qi and ty + qi, < q < 2ty for several choices 
of magnetic fields. 

We note here again that the above exponent seems 
to approach —1 for vanishing magnetic fields, while the 
correct exponent should be given by —1/2 at h — 0, cor- 
responding to the square-root singularity. As mentioned 
above, there is the crossover between these two behaviors. 
In this case, in the region (u — uj^q))^^ 1 ^/ 2 < l/r(2x), 
we can see a power-law behavior with the exponent 
X 2 = —1/2+ (2x — l)/2, while for higher energies, we 
can see the square root singularity. When h —> 0, this 
region with the exponent Xi is gradually shrunk, being 
replaced by the square-root behavior. 

FIG. 2. Momentum-dependent critical exponent Xi(q) 
and X%(q), which correspond to the lower edge of spectrum 
in ty — qL < q < ty + qL and < q < 7T — qi, ty + qL < q < 2ty. 
The magnetization is given by m — 0, 0.5 and 1 with the fixed 
U = 6. 
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C. Comparison with photoemission spectrum 



As shown above, the edge singularity in the dynamical 
charge correlation function for the ID Mott insulator is 
controlled by massless spinons subject to the phase shifts 
caused by the final-state interaction. We recall again 
that similar anomalous properties have been already re- 
ported in the study on the photoemission spectrum, for 
which one electron is removed from the system and the 
induced massless spinon excitations play a crucial role 
for the edge singularity. ]l6|— |l9[| Therefore, the origin of 
the anomalous critical behavior in both cases are essen- 
tially same in the sense that their class belongs to the 
Tomonaga-Luttingcr liquid with mobile impurities [ |34| . 
It should be noticed, however, that in the present case 
for the charge correlation function two mobile impurities 
are created as a particle-hole excitation, while for the 
photoemission spectra, only a single mobile impurity is 
created. So, the momentum-dependence of the critical 
exponents are quite different between these two cases. 



IV. SUMMARY 

We have studied the critical behavior in the dynamical 
charge correlation function for the ID Mott insulator. 
Anomalous critical exponents have been calculated for 
the edge singularity in the charge excitation spectrum for 
the half-filled Hubbard model, by combining the Bethe 
ansatz with finite-size scaling methods in CFT. It has 
been claimed that the final-state interaction between the 
charge and spin degrees of freedom plays an important 
role to produce the anomalous critical behavior depen- 
dent on the momentum. We have further pointed out 
that the universality class of the present system belongs 
to the Tomonaga-Luttinger liquid with two mobile impu- 
rities. 

In this paper, we have studied the charge density corre- 
lation function for the half-filled Hubbard model. Similar 
analysis can be straightforwardly applied to the dynam- 
ical spin correlation function for spin gapped systems 
in a metallic phase. In such cases, the roles played by 
the charge and spin degrees of freedom are interchanged; 
namely massless holons subject to the spinon phase shifts 
essentially determine the anomalous critical properties in 
the spin excitation spectrum. 
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